Delocalized Entanglement of Atoms in optical Lattices 
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We show how to detect and quantify entanglement of atoms in optical lattices in terms of cor- 
relations functions of the momentum distribution. These distributions can be measured directly in 
the experiments. We introduce two kinds of entanglement measures related to the position and the 
spin of the atoms. 
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Experiments on atoms in optical lattices have recently 
become very attractive playgrounds to investigate ba- 
sic issues in the context of quantum information theory 
[11, S S]- The high degree of control reached in those 
experiments should allow us to prepare a large variety 
of entangled states, to analyze their physical properties, 
and to verify and quantify the presence of entanglement. 
However, in order to carry out these investigations, a 
clear definition of entanglement with a clear physical 
meaning should be given, and ways to detect it should 
be explored Q ■ For atoms in optical lattices those ques- 
tions are far from being trivial, since one has to con- 
sider different degrees of freedom (such as atom numbers 
and internal levels) and take into account the presence of 
super-selection rules, as well as the fact that in most of 
those experiments we only have access to certain collec- 
tive properties. 

In this letter we will define and explore certain kinds of 
entanglement which are relevant in current experiments 
dealing with bosonic and fermionic atoms in optical lat- 
tices. We will concentrate on the entanglement proper- 
ties between different sites of the optical lattice (i.e. in 
second quantization), since in this case a clear meaning 
as a resource for quantum information tasks can be as- 
signed to those definitions. Furthermore, we will restrict 
ourselves here to the simplest case of bipartite entangle- 
ment. Since it is very hard in practice to address atoms 
at different lattice sites, we define the delocalized bipar- 
tite reduced density operators for a state p on an optical 
lattice by 



PAB 



m,m-{-x) 1 



(1) 



where P(m.m+x) denotes the restriction of p to the sites 
m and m -I- a; of the lattice. Since we do not want to rely 
on any form of addressability, we assume the lattice to 
be (approximately) infinite. To simplify matters, we con- 
sider here only a ID lattice, but everything holds as well 
when the state is defined on 2D or 3D lattices. We will 
study the entanglement of pab by means of experimental 
feasible collective measurements, that can be translated 
directly into expectation values of Pab- Note that, for 
this one to one correspondence between measurement re- 
sults and expectation values of paBi we define this state 



to be unnormalized. This definition of pab is useful for 
entanglement detection, however for a quantitative in- 
vestigation, we need to include the normalization in the 
definition. Due to the infinite number of lattice sites 
the trace of pab is infinite, such that a straight forward 
normalization fails to give any utilizable quantitative in- 
formation. Therefore we have to take into account that, 
for a finite number of atoms, most of the lattice sites are 
empty and do not contribute to any measurement. The 
obvious solution is to restrict the summation to a finite 
part of the lattice, such that pab can be normalized. 
But this makes only sense if we can ensure the atoms to 
be located in a relative small region of the lattice, what 
requires either some form of addressability or an extra 
assumption about the localization of the state p. A more 
general solution is to define the state as the (nor- 
malized) projection ol pab to the subspace where at least 
one atom is present, i.e., project out the zero atom sub- 
space of Pab- Even if the atoms are evenly spread over 
the whole lattice, as long as the number of atoms is small, 
the projection of pab has a finite trace. Unfortunately, 
we can not give the exact trace in terms of experimen- 
tal feasible measurements, but we can bound the trace 
which still allows us to derive lower bounds for the en- 
tanglement. Through the definition of pab and p'j^^ we 
are mixing different contributions of different sites, which 
implies that the entanglement we define will be somehow 
delocalized between different pairs of sites which are sep- 
arated by a distance x. Note, that pab depends on the 
distance x, which can be freely chosen. We will give lower 
bounds to the entanglement of formation for p'ab- Note, 
that due to the translationally symmetry of ([1]), we can 
never reach a maximally entangled state for p'ab- The 
maximally possible entanglement of formation that can 
be reached is 0.285. 7] 

We will consider two situations: firstly, that in which 
the atomic internal levels are not involved, and therefore 
we will deal with different occupation numbers; secondly, 
the one in which entanglement occurs between different 
internal states of the atoms in each site. As we will show, 
the first case is very simple to characterize and one can 
already claim that this kind of entanglement has been 
created in several experiments carried out so far. The 
second kind of entanglement is much subtler, and require 
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more sophisticated measurements in order to prove the 
existence of entanglement. 

The trapped atoms are described in second quantiza- 
tion by the annihilation and creation operators. We will 
assume two internal levels a and b for each atom, in which 
Qm, o-ln resp. bm, b\^ denote the annihilation and creation 
operators of atoms in level a resp. b at site m. 

Measurements: One kind of measurement that is fea- 
sible in experiments is to turn off the lattice potential 
and look at the density na{x,t) = ('0a(a;, i)^'(/'a(a;, i)) 
(resp. n},(x,t)) of the expanding atom cloud, where 
tpa{x,t) is the bosonic or fermionic field operator for 
the internal level a. Furthermore, one can mea- 
sure the density-density correlations Cabix,x',t) = 

(^'ipa{x,ty'ijja{x,t)'ipl{x',t)'il;b{x'^t)'^ and in a similar way 

Caa{x,x' ,t),Cbb{x,x' ,t) and Cbaix,x',t) i5-]. In a long 
time of flight approximation the new positions of the 
atoms can be detected and becomes proportional to the 
initial momentum distribution k w mx/{ht). Due to 
this relation we can measure the momentum distribu- 
tion via na{k) = limt^oo na{^^, t) and their correlations 
Cxy{k,k') = \\mt^ooCxy{^,^,t). This momentum 
distribution in second quantization is given by 



na(k) 



y^^Wn{k)wy, 



(2) 



and the density-density correlations by 

Cab{k,k')K (3) 

Wnik)Wm{k)*Wn>{k')Wm'{k')* ( 



n,m,n' ,m' 



where Wn{k) denotes the Fourier transformed Wannier 
function at time zero at site n. And in an analogous 
manner we get nb{k), Caa{k,k'),Cbb{k,k') and Cba{k,k'). 

Occupation number entanglement: We assume now the 
simple case, where we only have one type of atom, say 
that in level a. Due to the atom number conservation 
the possible product states in second quantization are 
restricted to be of the form 



, no, ni, n2. 



(4) 



because it is not allowed to have superpositions between 
states with different number of atoms. Here rii denotes 
the occupation number of sites i. We define a state to be 
entangled, if it can not decomposed into states of the form 
(j4|). Testing separability simplifies in this case to check 
whether a given state is diagonal in (j4|). Note, that en- 
tanglement defined with respect to such a super-selection 
rule is in general a less powerful resource for quantum in- 
formation tasks, because the entanglement properties can 
only be seen /used when having access to several copies 
of the state [6|] . To detect this kind of entanglement we 
look at the momentum distribution A simple obser- 
vation is that for all product states we get a 5{n = m) in 



the sum, because a\am maps any state of form ([4]) to a 
orthogonal one if n 7^ to. For this reason we get for any 
separable state that na{k) — N is just proportional to the 
total number operator and is independent from k, i.e., 
the momentum distribution is flat (up to the envelope 
Wannier- functions). Any non flat momentum distribu- 
tion indicates an entangled multipartite state, something 
that has been already observed in several experiments 
0. For a more quantitative statement about the entan- 
glement of ^ we look at the Fourier transformation of 
the momentum distribution of p 



:= Idk e-^^""^ {n. 



Here we have used, that multiplication by a phase e^^^'^ 
in momentum representation results in a shift in posi- 
tion, i.e., Wn{k)e~'^^^'^ = Wn+x{k) and that two at differ- 
ent places located Wannier functions are orthogonal, i.e., 
/ c?fc Wn+xW'!jy^ = S{n + x,m). x is taken to be an arbi- 
trarily integer and d denotes the lattice constant. For ^ 
we can give an interpretation in terms of an expectation 
value of the bipartite state pab ©: 



{Qx)p ■■= (4 



(6) 



Assume now the idealized case where the occupation 
number of every site is restricted to be either one or zero, 
defining this way exactly one qubit per site. Then pab is 
a two qubit density matrix and (^a\ag^ — (01|/5ab|10) 

is exactly one off-diagonal element, where |0) and |1) de- 
note empty or occupied sites. Due to the super-selection 
rules it is the only allowed off-diagonal element and 
defines the entanglement properties of the state. The 
state Pab can be decomposed into two parts. A separa- 
ble part spanned by the vectors |00) and the part 
spanned by the vectors |10), |01) that contains entangle- 
ment if (01|pab|10) 7^ 0. For this part we now want 
to estimate a lower bound for entanglement of forma- 
tion Q. Note, that we use here a definition for en- 
tanglement of formation respecting the super-selection 
rules. Assuming now a normalized state with given 
off-diagonal element A, it can easily be shown, that a 
pure state completely supported on the |10),|01) sub- 
space having the same off-diagonal element A gives us 
a lower bound to the entanglement of formation. Ex- 
ploiting this leads to a lower bound for the entanglement 

given by E^fip'AB) > ^[^[l- ^/T^W] ) , where 
S{x) = — a;log(a;) — (1 — x) log(l — x) denotes the von 
Neumann entropy. To estimate |A| for our we first 
have to find a bound for the trace of the unnormelized 
p'j^g . This can be given by 2 ( ) , since the the reduced 
densities sates in ([T|) cover two times the whole lattice. 

Therefore we can conclude that |A| > ^(n) ■ 

While the assumed restriction of maximal one atom 
per site matches perfectly for fermions, in the bosonic 
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case we can still give a bound if the following constraint 
can be guarantied, e.g., verified by further measurements 
0. The expected number of sites with more than one 
atom has to be smaller than e (N) and the maximally 
occupation number of one site has to be smaller than 
r. Under these conditions, the measurement result is 
still close to the off-diagonal element (01|pAs|10). The 
error coming from overpopulated sites can bounded by 
(2er + 4v^) (A'') (see Appendix) leading to 

(7) 




(2er + 4Vi) (A^) r)2 



{NY 



Internal level entanglement: We now consider the case, 
where we have two level atoms in the lattice. In the ideal 
situation we would have exactly one atom per site such 
that the internal levels realize one qubit. In this case pab 
is again a two qubit state without any restriction due 
to the conservation laws. To detect this stronger kind 
of entanglement it is not enough to look independently 
at the momentum distribution of level a and 6, but we 
have to look at the correlation 5] between momentum 
distributions Caa{k, k'), ci,b{k, k'), Cab{k, k') and Cba{k, k'). 
By properly chosen Fourier transformations in k and k' 
we define 



{Qf) 



dk dk' e 



ikxd , 



ik' xd 



E 



Cab{k,k') (8) 

(9) 



and in an analog way {Q°f') and {Q'^)- Here we again 
use the fact that the integrals over Wannier functions on 
different sites leads to delta functions. We furthermore 
assume a situation where we can restrict these sums to 
the case m — m! . 



{Qf) 



n h 

/ y \ "'m"'m+x^m+x^m 



(10) 



We will discuss later in the section dephasing how this 
condition can be realized by adding extra magnetic fields 
such that the m ^ m! terms vanish. 

The one atom per site case: To illustrate the basic 
idea, we assume now the idealized situation, where we 
have a state p for which we can ensure, that in every site 
is exactly one atom. Note, that in this case we can use 
()10|) without assuming any extra magnetic fields, because 
the terms with m ^ to' vanish already because of the 
assumption. Equation pop can now interpreted as the 
expectation value 



(11) 



of a bipartite density matrix pab as defined in ([T]). We 
want now to calculate the overlap of the state with 



a maximally entangled state, i.e., the fidelity /$(p^b) = 
where will be one of the Bell-state defined 

by 

0± = -^(|1O)a|O1)s±|O1)a|1O)b) 
V'± = -^(|10)a|10)b±|01)a|01)b) 

Here \10)a denotes the atom of Alice being in the a and 
|01)a being in the h level and in analog way for Bob. We 
claim that now, that 

(Qr + O, - ±((V.-),,, + (^+),,,) + (A^)/;2) 

where the (±) in distinguishes between the bosonic 
and the fermionic case. This can easily checked, by calcu- 
lating the expectation values for an arbitrary pure state 

1$) = Aoo|10)a|10)b + Aoi|10)a|01)s (13) 
+Aio|01)a|10)b + Aii|01)a|01)s. 

We get that (Qf )^ = AoiAto and (Q^")^ = X^^Xw such 
that the sum is equal to (— 10_)(</)_| -I- Fur- 
thermore we obtain that 



{QT) 



(14) 



\ ' PAB 



where the (±) distinguishes the bosonic from the 
fermionic case and in analogous manner for (Q^''). It 

is easily verified, that (o^a'^a'^'^b'^bI + \ ^\^a^\^b 



lAooP + lAiip - (V'-) + {'4>+) and (aVa + ^a^a) = 

{N)n 

From it is now easy to derive the fidelities /^^ for 
p'ab- 0^*3 to the one atom per site assumption we can 
bound the trace of the projected pab by {N) (instead of 
2 (TV)) leading to 



/0±(pTs)> 2 + 



2(iV) 



1 ±(Q°'' + Q''°)^ + (Q- + Q^''^A^), 
}^±{Pab) > 2 + 

for bosonic and fermionic case. Note, that we can get 
fidelities with further maximally entangled states by ap- 
plying a global unitary U ® ■ ■ ■ ®U to p. This translates 
to apply U ® U to pab- While the overlap with |(/)_) 
stays unchanged, we get that 

{4>+\U(^Upab{U(^U)^\4>+) = {4>u\pab\M, 

where \4)u) = 1 (8) {UU'^)\(j)+) is a maximally entangled 
state. In particular, we can get /^^ by properly chosen 
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U . The fidelities f^j^ and /^^ of a state directly gives us 
lower bounds to several entanglement measurements by 
comparing it to Bell-diagonal or Isotropic states e.g., 



Eos{p'ab) >S[\ 



1- Vl-(l-2/^, 



(15) 



The general case: We consider now the general case 
where the number of atoms per site is arbitrarily. In this 
case is it difficult to give a lower bound to {<t>±\pAB\4'±) i 
because the four Bell-states do not build a basis in the 
larger Hilbert-space and it is therefore not possible to 
derive the fidelity from (H^). But it still makes sense 
to define A -.^ ±{{(j}+) - (</>-)) + + in the 

larger Hilbcrt-space. Given a value of A it can be shown 
that the bound (fT5|) still holds with f^p^ now replaced by 
f<t>± ~ 1 ^■^^n though /^^ has now no direct interpre- 
tation as fidelity. Note, that A is a direct bound for the 
concurrence [1| and can be used itself to quantify the en- 
tanglement. Our goal is to give a lower bound for f^^ for 
a state p'^^ that we define in this case as the projection of 
PAB to the subspace with 2 or more atoms. We assume, 
that number of defective sites D should be bounded by 
D < e (N) . As defect counts every site having two or 
more atoms. Furthermore we assume the maximal occu- 
pation number of a and b to be less than r. A straight 
forward calculation (see Appendix) leads to 

(±(Q'''' + Q'"') + (2 - 4er2)iV - Q"" - Q^^) 



> 



2Wp 



for the bosonic and fermionic cases, which can be used 
to estimate a lower bound for p'^b- 

Dephasing with magnetic field: We discuss now the 
possibilities of eliminating the terms of ^ where m ^ m' 
by dephasing, i.e., by destroying any delocalization of 
atoms over several sites. If we write $ = J2k '^-R'I-^'^) in 
a product basis \K) = {kt^, . . . fc^, ktooi ■ ■ ■ where 
^(a/'') jjgjj^Qtes the number of atoms in site i in level a/6, 
then the unwanted terms that contribute to the sum are 
of the form 



(16) 



where m ^ m! . To give a nonzero value, \K) has to 
be equal to \K') although \K) has an extra atom in m 
and m' -\- X and a missing atom in m -|- cc and m', which 
implies defects in this pair of sites. We will use this dis- 
placed atoms and an additional inhomogeneous magnetic 
field to introduce some random phase to and there- 
with guarantee that these terms vanish. Let us assume 
a magnetic field, that is proportional to fc^, where k is 
the number of the lattice site. Then the state \K') gets, 
up to a global phase, a phase of e*('"+^) t+iva whereas 



\K') picks up e^'"'*+^(™'+^)'*. So \K) and \K') get a 
relative phase of e*^*^™ -m)xt By a, properly chosen 
set of times, we can randomize the phases such that for 
given X all terms with m ^ m! vanish in average, while 
the terms with m — m' stay unchanged. Note that even 
without magnetic field it seems to be quite unlikely, that 
all terms of the form (|16p sum up to a nonzero contri- 
bution. To give some non vanishing amount, the defects 
have to be correlated in a very unlikely way. In detail, 
one defect located at m or to + 2: has to be correlated 
to a defect located at to' or to' -|- x and in addition all 
these cases has to be correlated with each other. If we 
assume, that defects occur only randomly and therefore 
are uncorrelated, p6p already vanish. 

In conclusion, we have defined two figures of merit to 
quantify entanglement for atoms trapped in optical lat- 
tices, by only measuring density correlation functions of 
the expanded atomic cloud, without requiring any ad- 
dressability of the original lattice setup. One set of mea- 
surement data can be used to study entanglement at ar- 
bitrarily distances x. We discuss bounds in the cases 
where the defects in the lattice can be bounded. We ac- 
knowledges support from EU projects SCALA and DFG- 
Forschungsgruppe 635. 



[1] 
[2] 



[3] 



[4] 
[5] 

[6] 



[7] 



[9] 
[10] 



J. I. Cirac and P. Zoller, Phys. Today 57, 38 (2004) 

D. Jaksch, H.-J. Briegel, J. I. Cirac, C. W. Gardiner, 
and P. Zoller, Phys. Rev. Lett. 82, 1975 (1999); G. K. 
Brennen, I. H. Deutsch, and C. J. Williams, Phys. Rev. 
A 65, 022313 (2002); J. Mompart, K. Eckert, W. Ertmer, 
G. Birkl, and M. Lewenstein, Phys. Rev. Lett. 90, 147901 
(2003); K. G. H. VoUbrecht, E. Solano, and J. I. Cirac, 
Phys. Rev. Lett. 93, 220502 (2004) 

O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. 
Hiinsch, and I. Bloch, Nature, 425, 937-940 (2003); 
M. Greiner, I. Bloch, O. Mandel,T. W. Hansch and T. 
Esslinger, Phys. Rev. Lett. 87, 160405 (2001); H. Moritz, 
T. Stoferle, M. Kohl and T. Esslinger, Phys. Rev. Lett. 
91, 250402 (2003); M. Kohl, T. Stoferle, H. Moritz, C. 
Schori, T. Esslinger, Applied Physics B 79, 1009 (2004). 

G. Toth and O. GiLhne, Phys. B 82, 237 (2006). 

E. Altman, E. Demler, and M. Lukin, Phys. Rev. A 70, 
013603 (2004). 

Norbert Schuch, Frank Verstraete, and J. Ignacio 
Cirac,Phys. Rev. A 70, 042310 (2004); S. D. Bartlett, 

H. M. Wiseman, Phys. Rev. Lett. 91, 097903 (2003); 

U. V. Poulsen, T. Meyer, D. Bruss, and M. Lewenstein 
Phys. Rev. A 73, 052326 (2006); K. M. O' Connor and 
W. K. Wootters, Phys. Rev. A, 63, 052302 (2001). 
C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, W. K. 
Wootters, Phys. Rev. A54 (1996) 3824-3851; K.G. H. 
VoUbrecht, R.F. Werner, Phys. Rev. A 64, 062307(2001); 
W. K. Wootters, Phys. Rev. Lett. 80 (1998) 2245-2248. 
P. Rabl, A. J. Daley, P. O. Fedichev, J. I. Cirac, and P. 
Zoller, Phys. Rev. Lett. 91, 110403 (2003). 
We assumed, that atoms in level a and b gets the same 
phase independent of the internal state. 



5 



Appendix 

Calculating errors for the off-diagonal element: pab 
can be written as a direct sum X^nP"' where p„ is the 
n-atom subspace. The off-diagonal element (01|pyis|10) 
is exactly given by {Qx)p-^- So we have to bound all the 
absolute values of {Qx)p for n > 1. Note, that all p„ for 
n > 2 has at least one defect. Therefore we can bound 
the trace by two times the number of defects (2e (N)) and 
therefore the total contribution by 2e (N) r. We need to 
take twice the number of defects, because the reduced 
states in ([T]) cover two times the whole lattice. More 
complicated is to find a bound p2, because 2 atoms do 
not imply automatically a defect. We therefore look at a 
pure state supported on p2 

|^2) = Ao|ll) + Ai|20)+A2|02). (17) 

(Qa;)^2 easily calculates to (AqA^ + AqA2)V2, with 
J2i |AiP = 1. We now try to find an upper bound for 
its absolute value. First lets assume, that all Xi are real 
and positive leading to Ao(Ai + A2)V2- For given Ao this 
is maximized for Ai = A2. Defining A' = "^vAf + Al, 

we get (Q^)^^ < Ao2A'V2 < 2A'V2 < 2V'Af+A| as up- 
per bound. Now, we use that ^/x < -^^x + for any 



positive parameter R, which gives us {Qx)^^ < "^(A? + 
A2) + \/R. Here (Af + Aj) is now exactly the probability 
for an defect. For the unnormclized p2 = Tlii1i\'^2){'^2\ 
we get {Qx)p^ ^Vn^^ + ^2)* + Vi? trip2). The 
sum J2i <Zi(Ai -I- A|)j is now the probabihty to found a de- 
fect in p2, that can be bounded by 2e {N). Furthermore, 
the trace of p2 can be bounded by tr{p2) < {N), leading 
to {Qx)p^ < (^^^ ^ V^)2 (N). Now we a free to choose 
R = y/e and get finally {Qx)p, < 4^{N). So we get 
|(01|pab|10)| = I {Qx)p, I > I {Qx)p I - (4V^ + 2er) (N) 

Calculating errors for fidelity: We want to bound A = 
2/0± — 1 = ± {4>+4>+) + (V'+?A+). Given the state pab = 
Pn written as direct sum. We are now interested in 
^PAB which is equal to Ap^, because A is supported on 
the 2-atom subspace. Since holds also on the bigger 
2-atom subspace 2/0^(^2) — 1 equals Ap^. Ap^ gives no 
contribution and can therefore neglected. For this reason, 
we only have to take defects with two or more atoms 
per site into account. We therefore have only to bound 
2/0± [p'^) — Ij where p-^ = X]n>2 Note, that this sum 
runs by assumption only up to r^. The trace of p^ is 
smaller than 2D < 2e {N) , since for n > 2 there exists 
always a defect and 2/^^ — 1 can be bounded by 4r^ for 
bosons and by 4 for fermions. 



